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I.  INTRODUCTION 
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gaOEPlIP  PA®  Ei  |HT»MOT  F1U® 


This  work  is  a  continuation  of  the  efforts  of  Mullinix  and  Smith  [  I  land  Vandiver  etal. 
1 2].  Their  work  involved  the  determination  of  stress  intensity  factors  for  homogeneous 
cylinders  loaded  under  one  of  three  conditions:  internal  pressure,  bending,  extension  and 
having  part  circular  longitudinal  or  transverse  cracks  in  the  outer  wall.  The  present  work 
considers  the  problem  of  finding  stress  intensity  factors  under  the  combined  loading  action  of 
internal  pressurization  and  extension.  Transverse  flaws  in  cylinders  were  considered  for  this 
effort. 


Although  flaws  on  the  surface  of  cylinders  usually  occur  as  semi-elliptical  in  shape,  the 
flaws  in  this  work  had  to  be  made  mechanically.  This  was  accomplished  by  cutting  part- 
circular  flaws  with  circular  saw  blades  to  simulate  an  elliptical  crack.  The  approximation  has 
been  made  before  and  does  not  appear  to  offer  any  serious  error.  In  this  effort  the  assumed 
stress  field  at  the  crack  tip  border  is  defined  as  a  50-50  stress  mixture  derived  from  internal 
cylinder  pressurization  and  extensional  loading.  Fifty  percent  of  the  applied  stress  was 
obtained  from  internal  pressurization  while  the  remaining  fifty  percent  was  obtained  from 
extension. 

U.  THEORY 


The  geometry  for  the  part-circular  transverse  flaw  may  be  described  by  the  intersection 
of  a  circular  element  representing  the  flaw  boundary  with  a  hollow  cylinder.  Figures  1  and  2 
illustrate  this  geometry.  For  the  opening  mode  of  deformation,  the  stress  distribution  near  the 
part-circular  crack  and  in  a  plane  perpendicular  to  the  crackfront  is  given  as  [3,4,5] 
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ln  Equation  (1),  Ki  is  the  stress  intensity  factor  and  the  coordinates  for  the  crack  border 
are  shown  in  Figure  3.  The  effect  on  the  stress  field  of  the  crack  border  curvature  as  well  as  the 
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location  and  shape  of  other  boundaries  is  reflected  in  the  magnitude  of  Ki,  the  stress  intensity 
factor.  It  is  assumed  that  the  stresses  singular  in  r  are  much  larger  than  terms  regular  in  r  very 
near  the  crack  tip.  Since  experimental  measurements  must  be  made  away  from  the  crack  tip 
where  regular  terms  may  make  a  significant  contribution  to  the  total  stress  field  then  another 
way  of  determining  Ki  must  be  found.  Irwin  J 6 J  developed  an  approximation  for  the  regular 
terms  by  assuming  a  uniform  stress  field.  o„n,  was  superimposed  at  the  crack  tip  and  parallel  to 
the  crack  plane.  With  his  approximation  the  local  stress  field  is  given  as 
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Oon  does  not  affect  the  singular  stress  field  but  does  alter  the  isochromatic  fringe  pattern  which 
is  proportional  to  the  maximum  in-plane  shearing  stress.  The  maximum  shearing  stress,  rma„is 
usually  determined  readily  from  photoelasticity. 

From  Irwin’s  stress  equations,  the  maximum  shearing  stress  in  the  plane  perpendicular 
to  the  crack  front,  y-n,  can  be  obtained  from 
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In  photoelasticity,  the  maximum  shearing  stress  from  the  stress  optic  law  is 


T 

max 


£N 

2t 


(5) 


where, 

t  =  thickness  of  the  specimen  measured  parallel  to  the  direction  of  light  propagation. 
N  =  isochromatic  fringe  order, 
f  =  photoelastic  fringe  constant  of  the  material. 


In  practice,  Tmslx,  is  measured  along  the  line  i/»  =  jr/2  where  the  maximum  shearing  stress  is 
known  to  be  large.  Simplification  of  Equation  (4)  with  «/»  =  7r / 2  results  in  the  equation 
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Only  the  stress  intensity  factor  Ki  is  used  in  fracture  criteria.  Smith  et  al.  |7]  found  that  by 
solving  Equation  (6)  for  rma»  and  truncating  the  results  to  the  same  order  in  r  as  Equation  (2) 
that 


T 

max 

If  an  apparent  value  of  K,  Kap  is  defined  as, 

K  =  \|  8irr  x  ,  (8) 

ap  N  max  '  1  ' 

then  Equation  (7)  can  be  written  as 

Kap  =  h  +  x/®?r  °on  (9) 

Equation  (9)  shows  that  in  the  region  dominated  by  the  singular  stresses  that  there  is  a  linear 
relationship  between  the  apparent  K  and  the  square  root  of  r.  In  determining  Ki.  the  values  or 
Kap  are  plotted  versus  r1 :  for  a  photoelastic  slice  specimen.  Data  points  which  fall  on  a  straight 
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line  are  selected  while  all  others  are  rejected.  A  least-squares  straight-line  curve  fit  is  then  given 
to  the  selected  points  and  the  value  of  Ki  is  determined  by  taking  the  value  of  K,pat  r  =  0,  since 
Kap  =  Ki  at  r  ~  0.  Figure  4  gives  an  example  for  a  typical  photoelastic  slice  specimen. 


For  clarity,  the  results  of  this  experimental  effort  are  compared  with  those  of  Reference 
[2|.  To  compute  the  overall  stress  level  for  experimentally  subjecting  a  cylinder.  Reference  [5] 
was  consulted.  Thresher  and  Smith  generated  graphs  of  stress  intensity  factors  for  surface 
cracks  in  finite  solids.  Using  their  information  along  with  the  maximum  allowable  working 
stress  in  the  photoelastic  material  a  determination  of  om,  the  nominal  cylinder-wall  stress,  was 
made.  From,  am,  the  maximum  working  internal  pressure  P,  of  the  cylinder  was  obtained  from 


where 


(10) 


T  =  Wall  thickness  of  the  cylinder. 

Rc  =  Radius  of  cylinder  measured  to  the  center  of  the  cylinder  wall. 

For  the  case  of  pure  internal  pressure  loading  of  a  cylinder  reported  in  [2], 


P. 

l 


(11) 


To  compute  the  extensional  loading  P  for  the  cylinder  with  a  50-50  stress  mixture, 


P  =  1  0  A 
2  ^  ^ 


(12) 


while  for  the  case  reported  for  pure  extensional  loading  in  [2J, 
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a 
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A 
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(13) 


All  stress  levels  were  comparable  since  their  ultimate  determination  was  from  the  same  source. 
Thresher  and  Smith  [5]. 

III.  EXPERIMENTATION 


The  determination  of  stress  intensity  factors  for  cylinders  loaded  in  internal  pressure 
and  extension  followed  the  work  given  in  References  1 1  ]  and  |2|  which  used  three  dimensional 
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photoelasticity.  A  series  of  seven  combined  loading  tests  were  performed.  The  photoelastic 
material,  Hysol  CP5-4290,  was  cast  by  the  Hysol  Corporation,  Olean,  New  York,  and  used  in 
the  experimentation.  T  he  cylinders  are  nominally  5.875  inch  in  outside  diameter  with  a  0.75- 
inch  wall  thickness.  All  the  specimens  had  flaws  oriented  transverse  to  the  cylinder  axis.  These 
flaws  were  machined  with  a  circular  saw  blade  0.006-inch  thick.  Blade  radii  of  0.875  and  1 .500 
inch  were  used  to  produce  flaws  of  two  different  sizes. 

For  the  seven  tests,  the  internal  pressure  and  extension  loads  were  determined  using 
Equations  (10)  and  (12).  Figure  5  illustrates  the  apparatus  for  generating  the  internal  gas 
pressure  and  extensional  loading.  The  uniform  tension  load  was  supplied  by  hanging  dead 
weights  on  the  cylinder.  The  internal  pressure  load  was  supplied  by  compressed  air  passing 
through  a  regulator.  The  gas  pressure  was  measured  by  a  Mercury  manometer. 

After  the  surface  flaw  was  machined,  the  cylinder  was  annealed  by  thermal  soaking  at 
280  degrees  Fahrenheit  for  six  hours  followed  by  cooling  at  the  rate  of  one  degree  Fahrenheit 
per  hour.  The  stress  freezing  of  the  models  was  accomplished  using  the  same  thermal  cycle  as 
for  annealing  except  under  loading  conditions.  After  the  stress  freezing  cycle,  slices 
perpendicular  to  the  crack  border  were  removed  from  the  model  by  means  of  band  saw.  Figure 
6  illustrates  a  few  of  the  various  angles  at  which  slices  were  taken.  The  number  of  slices  varied 
from  test  to  test  depending  on  the  flaw  size.  Each  slice  was  polished  with  sandpaper.  The  CP5- 
4290  material’s  fringe  constant  was  obtained  from  previous  beam  tests. 

To  improve  resolution  for  analysis,  slices  were  placed  in  an  oil  bath  consisting  of  75.5 
percent  by  volume  of  Halowax  oil  and  24.5  percent  mineral  oil.  Since  the  indices  of  refraction 
of  the  oil  and  CP5-4290  were  the  same,  light  scatter  was  minimized.  The  slices  were  observed  in 
a  comparator  polariscope  at  10X  magnification. 

By  means  of  an  XY-table  on  the  polariscope.  points  on  the  slices  could  be  located  to 
within  ±0.0001  inch.  Fractional  fringe  orders  were  obtained  using  Tardy  compensation. 

IV.  RESULTS  AND  DISCUSSION 

The  specimen  test  parameters  and  dimensions  for  the  combined  uniaxial  tension  and 
internal  pressure  loading  tests  are  indicated  in  Table  3.  Tables  I  and  2  were  reproduced  from 
Reference  [ 2 ]  for  comparison  purposes  and  are  for  the  separate  loading  cases.  Figures  7-10 are 
graphs  of  the  non-dimensional  stress  intensity  factor  versus  slice  angle  for  transverse  flawed 
cylinders  loaded  either  in  uniaxial  tension  or  internal  pressure.  Figures  11-17  axe  graphs  of  the 
non-dimensional  stress  intensity  factor  versus  slice  angle  for  transverse  flawed  cylinders 
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loaded  in  combined  uniaxial  tension  and  internal  pressure.  The  combined  loading  case  data 
compares  favorably  with  the  separate  uniaxial  and  internal  pressure  load  cases.  In  general,  the 
data  follows  the  same  experimental  trends  and  the  stress  intensity  factors  obtained  for  the 
combined  loading  case  falls  within  the  range  of  the  individual  load  cases. 

V.  SUMMARY  AND  CONCLUSIONS 

A  series  of  seven  separate  tests  were  conducted  in  which  part  circular  flaws  simulating 
natural  elliptical  cracks  were  machined  into  Hysol  CP5-4290  cylinders.  The  cylinders  were 
subjected  to  combined  internal  pressure  and  uniaxial  tension  loading.  A  photoelastic  stress 
freezing  cycle  was  conducted  for  each  cylinder.  Following  the  stress  freezing  cycle  each 
cylinder  was  sliced  and  analyzed  using  a  polariscope  with  Tardy  compensation.  The  stress 
intensity  factors  were  shown  plotted  versus  slice  angle  and  were  compared  with  a  previous  set 
of  tests  reported  in  Reference  [2]  The  results  indicate  that  the  more  complicated  combined 
loading  case  produces  results  comparable  to  the  separate  loading  cases.  It  appears  that  linear 
superposition  of  solutions  for  each  stress  intensity  factor  case  (i.e.,  uniaxial  loading  or  internal 
pressure  loading)  is  valid. 


p 

Figure  1.  Transverse  flaw  loading  geometry  for  a  hollow  cylinder 


Figure  2.  Notation  for  the  part-circular  surface  flaw. 


Figure  4.  Typical  set  of  slice  data,  illustrating  the  determination  of  K,. 
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Figure  5.  Schematic  configuration  for  internal  pressure  and/or  extension  loading  of 
cylinders. 
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Figure  8.  Stress  Intensity  factor  versus  0  for  the  transversely-flawed  cylinder 
loaded  in  uniaxial  tension  (R  =  1.500  inch). 


17 


CIRCUMFERENTIAL  FLAW  WITH  INTERNAL  PRESSURE 

R  =  0.875  INCH 


Figure  9.  Stress  intensity  factor  versus  0  for  the  transversely-flawed  cylinder 
loaded  with  internal  pressure  (R  =  0.875  inch). 


(JifA)  uio(*/3) 


CIRCUMFERENTIAL  FLAW  WITH  INTERNAL  PRESSURE 
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Figure  10.  Stress  intensity  factor  versus  0  for  the  transversely-flawed  cylinder 
loaded  with  internal  pressure  (R  =  1.500  inch). 
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Figure  1 1 .  Stress  Intensity  factor  versus  0  for  the  transversely-flawed  cylinder 

loaded  with  Internal  pressure  and  extension (R  =  0.875  Inch,  A/T  =  .200). 
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Figure  12. 


TRANSVERSE  FLAW  WITH  COMBINED  LOADING 
(R  =  0.875  INCH) 
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Stress  intensity  factor  versus  6  for  the  transversely- flawed. cylinder 
loaded  with  Internal  pressure  and  extension  (R  =  0.875  inch,  A/T  =  .400). 
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Figure  13.  Stress  intensity  (actor  versus  8  (or  the  transversely-tlawed  cylinder 

loaded  with  internal  pressure  and  extension  (R  =  0.875  inch,  A/T  =  .500). 
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Stress  intensity  (actor  versus  8  for  the  transversely-flawed  cylinder 
loaded  with  internal  pressure  and  extension  (R  =  0.875  inch,  A/T  =  .600). 
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Stress  intensity  factor  versus  6  for  the  transversely-flawed  cylinder 
loaded  with  Internal  pressure  and  extension  (R  =  1.500,  inch,  A/T  =  .200). 


24 


(2Ar)om  (V»R) 


1.5 


TRANSVERSE  FLAW  WITH  COMBINED  LOADING 
(R  =  1.500  INCH) 


1.01 


& 


0-5H 


0.0+ 

0 


Figure  16. 


—r 

10 


~r 

15 


e(DEG) 


T 

20 


Stress  intensity  factor  versus  6  for  the  transversely-ftawed  cy Under 
loaded  with  Internal  pressure  and  extension  (R  =  1.500  Inch,  A/T  =  .400). 
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Figure  17.  Stress  intensity  factor  versus  9  for  the  transversely-flawed  cylinder 

loaded  with  Internal  pressure  and  extension  (R  =  1.500  Inch,  A/T  =  .600). 
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I  his  appendix  contains  a  summary  oi  all  the  test  data  presented  in  this  report.  The  test 
nomenclature  is  as  follows: 


T  s  Extension  loading  test. 

PC  =  Internal  pressure  loading  test. 

=  Combined  internal  pressure  and  extension  loading  test. 


PE 


The  test  nomenclature  follows  the  test  number.  Tor  example.  I  PE  refers  to  the  first  test  run  of  a 
flaw  in  the  transverse  orientation  subjected  to  combined  internal  pressure  and  extension 
loading. 


TEST  NO.  IT 


A/T  =  0.206 

R  --  0.675  Inch 

°m  *  10.74  pal 

SLICE  NO. 

ANGLE 

(DEG) 

*1 

Kl 

2/«  ffm<v£R) 

1 

0 

7  1788 

0  633 

2 

10 

6  0309 

0  532 

3 

15 

46534 

0410 

4 

20 

4  4673 

0  394 

TEST  NO.  2T 


A/T  =  0.390 

R  -  0.875  Inch 

"m  =  10.32  pal 

Kl 

ANGLE 

2/»  <rm(  v»R) 

SLICE  NO 

(DEG) 

Kl 

1 

0 

10  7465 

0  986 

2 

10 

10  6959 

0  981 

3 

15 

9  9365 

0  912 

4 

20 

9  3772 

0  860 

5 

25 

8.7716 

0  805 

I 


TEST  NO.  3T 

A/T  -  0.551 

R  0.875  Inch 

°m==9- 

Kl 

ANGLE 

SLICE  NO 

(DEG) 

K, 

2/»  ffm(v£R) 

' 

1 

0 

11  9434 

1  164 

2 

15 

12  6552 

1  233 

3 

20 

11.3063 

1.102 

4 

25 

12  4901 

1  217 

L 

5 

30 

12  2604 

1  195 

TEST  NO.  4T 

A/T  =  0.674 

R  -  0.875  Inch 

3BA 

r 

Kl 

ANGLE 

SLICE  NO 

(DEG) 

K, 

2/» 

TEST  NO.  2PC 
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APPENDIX  B 
COMPUTER  CODE 
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I  he  computer  code  shown  on  the  following  pages  was  used  to  analyze  the  photoelastic 
data.  I  he  code  uses  the  technique  presented  in  Section  II.  With  at  least  squares  straight-line 
curse  lit  of  the  experimental  data. 


C - PHOIOEI  ASTICITY  CODE-CYLINDERS 

DIMENSION  AN(50),AR(50).AT(50).AK(50),ID(50) 
WRIIEI5.22) 

22  FORMATE  NO.  OF  SLICES?') 

READ(5,I)  N 
1  FORMATf  13) 

DO  19  1=1, N.  I 
M  =  I4 
F=L56 
WRITE(5,25) 

25  FORMAT!'  SLICE  THICKNESS?') 

READ(5,26)  T 

26  FORM  AT(  F  10.0) 

WR1TE(5,27) 

27  FORMAT!'  INPUT  N-F6.0') 

DO  4  J=LM, I 
READ(5,3)  AN(J) 

3  FORM  AT(F6.0) 

IF(J.LE.9)  AR(J)=FLOAT!J-l)*.005+.0l0 
IF(J.GT.9)  AR(J)=FLOAT(J-9)*.0l0+.050 

4  CONTINUE 
AMAX=0. 

DO  5  J=LM,1 
AT(J)=F*AN(J)/(2.*T) 

AK(J)=AT(J)*SQRT(8.*3. 14159*AR(J)) 

IF( AK(J).GT.AM AX)  AM AX=AK(J) 

5  AR(J)=SQRT(  AR(J)) 

6  CALL  IPOKE("I704IO,"1) 

CALL  IPOKE("I70410,"0) 

DO  7  J=I,M,1 

IY=1NT(AK(J)*I000./ AM  AX) 

IX=1NT(AR(J)*  1000./  A  R(M)) 

CALL  IPOKE("I704I2,IX) 

CALL  1POKE("I704I4,IY) 

7  CALL  IPOKE(''I70414,"0) 

ITEST=!PEEK("1 77570) 

IF(ITEST.EQ.O)  GO  LO  6 
WRITE(5.8) 


s 


pi 


9 


10 
1 1 

20 


13 

21 


14 


15 


16 

18 

17 

19 


FORMATE  NO.  ()!  DELETE  SAMPLES?  ) 

RLAIX5.9)  NO 

FORM  ATM3) 

lE(ND.EQ.O)  GOTO  20 

DO  1 1  J=I,ND.I 

READ(5. 10)  1D(.I) 

FORMAT(13) 

CONTINUE 

CONTINUE 

Xl=0. 

X2=0. 

Y0=0. 

Y  1=0. 

DO  14  J=I,M,I 
IF(ND.EQ.O)  GOTO  21 
DO  13  K=I,ND,I 
1F(J.EQTD(K))  GOTO  14 
XI=XI+AR(J) 

X2=X2+AR(J  )**2. 

Y0=Y0+AK(J) 

Y1=Y1+AK(J)*AR(J) 

CONTINUE 

AK  l=(X2*YO-Yl*XI  )/(FLOAT(M-ND)*X2-Xl*Xl) 
WRITE(5,I5)  I 

FORMATS  SLICE  NO.  =',I3) 

DO  18  J=1,M,1 
AR(J)=AR(J)**2. 

WRITE  (5,16)  AR(J),AN(J),AT(J),AK(J) 

FORMATE  R=',F10.4,5X,'N='F10.4,5X, 

ITM  AX=',F10.4,5X,'KAPP='FI0.4) 

CONTINUE 
WRITE(5, 17)  AK  1 
FORM  AT('  K1=',F14.4) 

CONTINUE 

STOP 

END 


I 


44 


LIST  OF  SYMBOLS 


A 

A, 

a_ 

D 

D 

d 

f 

K.i,Sn 

N 

P 

P, 

R 
Rt 
r,  <l> 

T 

t 

y,  n,  0 

O  m 
Oon 

°V»  On,  Oq 
Tny,  Tnff,  TyQ 

T  max 
0max 


Crack  depth  at  deepest  point  (semi-minor  diameter  for  elliptical  flaw) 

Nominal  cross-sectional  area  of  cylinder 

Half-length  of  crack  on  outside  surface  of  wall 

Distance  from  center  of  circular  flaw  to  surface  of  the  wall 

Ratio  of  Distance  D  to  radius  of  circular  flaw  R 

Distance  from  center  of  cylinder  to  center  of  circular  flaw 

Photoelastic  fringe  constant 

Mode  1  stress  intensity  factor 

Isochromatic  fringe  order 

Total  load  on  cylinder  loaded  in  tension 

Internal  cylinder  pressure 

Radius  of  circular  flaw 

Radius  of  cylinder  measured  to  center  of  cylinder  wall 

Polar  coordinates  centered  at  crack  tip 

Wall  thickness  of  the  cylinder 

Thickness  of  a  slice  analyzed 

Coordinate  system  shown  in  Figure  3 

Nominal  cylinder-wall  stress 

Uniform  stress  at  the  crack  tip 

Normal  stress  components 

Shear  stress  components 

Maximum  shearing  stress  in  the  plane  perpendicular  to  the  crack  border 
Maximum  flaw  angle  (Figure  2) 
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